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^ . 1 Abstract 

Charged- current single pion production in scattering off ^^C is investigated for neutrino 
energies up to 1 GeV. An impact of nuclear effects with in-medium modifications of the 
. A(1232) resonance properties as well as an effective field theory nonresonant background 
i contribution are discussed. Dependence of the fraction of A(1232) decays into npnh states 
^ \ on incident neutrino energy is estimated. A model of Nieves et al. [TJ is further developed 
\Q • by performing exact integration avoiding several approximations. The effect of exact inte- 
^ ■ gration is investigated both for double-differential and total neutrino-nucleus cross sections. 
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2 Introduction 



^ ■ There has been a lot of effort to understand better the single pion production (SPP) reactions 
■ in neutrino-nucleon and neutrino-nucleus scattering. Motivations for these studies come 
from neutrino oscillation experiments and their demand to reduce systematic errors. In a 
few GeV energy region characteristic for experiments like T2K, MINOS, NOvA, MiniBooNE 
and MicroBooNE the SPP channels account for a large fraction of the cross section (at 1 GeV 
on an isoscalar target ~ 36%). 

In the neutrino experiments one typically measures charged current quasi-elastic (CCQE) 
events. The signal events have no pions in the final state. In the case of SPP reaction on nu- 
clei there is a significant probability that a pion produced on a bound nucleon gets absorbed, 
and such events contribute to the CCQE background. The MiniBooNE experiment has put 
a lot of effort to develop methods to estimate and subtract that background from the CCQE 
data sample [2J. A data/MC correction function was introduced based on the investigation 
of events with a single pion in a final state. However, it is not clear how legitimate was to 
assume that the same function can be applied to hypothetical pion absorption events. A 
necessity to use the large correction function can be an indication that the pion production 
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model implemented in NUANCE Monte Carlo event generator suffers from deficiencies [3]. 
Another well known instance of relevance of pion production channels is neutral current 7r° 
production. Neutral pions give rise to events which mimic — )■ z/g signal (it happens if one 
of the two photons from the 7r° decay remains unidentified). In the nuclear medium interme- 
diate virtual resonance states, leading normally to single pion production, may get absorbed 
by nucleons. The pionless A decays events contribute to multinucleon ejection final states 
and can be confused with the genuine CCQE events because the knocked out protons are 
usually not analyzed at all. The multinucleon knock-out contamination has an impact on 
the neutrino energy reconstruction and should be accounted for in the neutrino oscillation 
experimental analysis [4J. Weak single pion production processes are also important for the 
hadronic physics. They provide a valuable information both on the dynamical structure of 
the nucleon resonances and nonresonant contributions. This information is complementary 
to what is known from the electro- and phot opro duct ion studies. This topic was studied in 
the MAINZ, BONN and TJNAF laboratories. The results give a good insight on the elec- 
tromagnetic structure of nucleon resonances production, see for example the latest paper on 
the subject ^ and the underlying unitary isobar model for pion electroproduction [6j. The 
above mentioned analysis includes also a variety of low lying resonances beyond A(1232). 
Three of them: Pii(1440), Di3(1525) and 5*11(1535) may be relevant for the understanding 
of pion production process in the neutrino experiments, like T2K or NOvA. Electro- and 
photoproduction experiments can also serve as a valuable source of information on final state 
interactions (FSI) effects which are universal for all pion production experiments though a 
number of such studies is limited. 

Recent experimental results on the charge current SPP reactions come mainly from the 
K2K ([7], [8]) and MiniBooNE experiments ([9], [in])- Unfortunately, the analysis of the 
underlying fundamental physical processes of pion production on nucleons is obscured by 
nuclear effects. There is an important impact of the nuclear medium on a primary interaction 
as well as on a redistribution of physical channels by FSI inside the target nucleus. The 
nuclear physics uncertainties are so large that MiniBooNE did not attempt to measure 
the characteristics of neutrino- nnc/eon SPP process and published the cross sections results 
with all the nuclear effects included (the signal events are those with a single pion leaving 
a nucleus). A more valid information about the nucleon- A transition form factors and 
the effects of the nonresonant background can be extracted from the data obtained on light 
targets like hydrogen or deuterium. The results from two old low statistics dedicated bubble 
chamber experiments, ANL ([11], [12]) and BNL ([13]), are still used in the verification of 
dynamical pion production models and weak N — A transition form-factors (see for example 
U4J, ll5j). 

In the theoretical discussions of weak SPP two main concerns are: description of the 
non-resonant background and handling of the nuclear effects. For the neutrino energies of 
about 1 GeV a dominant pion production channel is that via A excitation. However, as 
clearly follows from the ANL and BNL experiments, there is an important non-resonant 
contribution that cannot be neglected though much harder to include in the computations. 
The nuclear effects include the A in-medium self energy which leads to a substantial fraction 
of pionless decays. An impact of the nuclear medium on the vector and axial A excitation 
form-factors is unknown and currently there is no method of measuring it. The problem of 
charged- current SPP on nuclei assuming A dominance model with many-body effects from 
[16j has been addressed in [IT]. The computations have shown a significant reduction of 
the pion production cross-section due to the in-medium effects. The fraction of A pionless 
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decays has a rather mild dependence on the incident neutrino energy. All these calculations 
did not include any kind of nonresonant background. An assumption of the constant fraction 
of pionless A decay is implemented in neutrino Monte Carlo event generators (NUANCE, 
NEUT) and the nonresonant dynamics is modelled in a simplified manner within an old 
fashioned Rein-Sehgal model which is known to suffer from other deficiences as well [3j. In 
this paper the full model of weak SPP on nuclei based on [14J is used for the neutrino- 
nucleus scattering following the approach of [1]. The impact of the nonresonant background 
is discussed. A lower bound for np — nh contribution coming from the pionless A decays 
is estimated. The goal of this paper is to present the predictions from the sophisticated 
theoretical model for SPP in such away that they can be used in the evaluation of the 
systematic errors by experimental groups. An impact of various ingredients of the model 
on the final results will be presented as well. 

Our model is based on the papers by Nieves et al.([l3], [1]) and we aim to further develop 
their approach. Our most important contribution is a prescription how to perform many 
integrals in an exact way. Thanks to that we avoid not easy to control approximations. 
It turns out that the approximations used in [1] do not work well in the case of double- 
differential cross sections, while for total cross sections they produce results close to the 
exact ones. We have found also, that the assumption of a constant fraction of pionless A 
decay cannot be applied for experiments with large flux contribution from Ey < 1 [GeV] 
and that the ratio of muon to electron (anti-)neutrino total crosss section does not depend 
on the medium modifications of A(1232) resonance. 

The paper is organized as follows: in Section [3] we discuss the general formalism of 
SPP on atomic nuclei. The dynamical model of SPP is reviewed in Section 13.11 and the 
nuclear medium effects are discussed in Section 13.1.11 In Section H] we briefly introduce the 
numerical procedures and in Section [S] we present our main results that are then discussed 
in Section E] 



3 Theoretical Description of Pion Neutrinoproduction 
on Atomic Nuclei 

The theoretical approach presented in this paper is based on the general scheme described in 
[1]. The basic cross-section formula for the electromagnetic or weak charged-current lepton 
inclusive differential cross section is: 



dn'dE' 



Fi(Q")j^ I d'vL.yW^'^ipir)) (3.1) 
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G|. cos2 Oc 



electrons 



; (anti)neutrinos 

+ ^'fi^iy ~ 9tiM'-, electrons 
^ij-^'u + Uu - 9i^M' ± ^(^^luaf}l"^l^ , (anti)neutrinos 



For the weak interactions the Fermi contact is Gp = 1.1664* lO^^^/Mel^^ and the cosine of 
Cabbibo angle is cos(Oc7) = 0.974. Furthermore and /'^ denote initial/final lepton four- 
momenta, = —Q^ = (/ — /')^ is the squared four- momentum transfer. In the laboratory 
frame we assume the momentum transfer to be directed along the Z-axis and the scatttering 
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to take place in the X-Z plane. The local density approximation is adopted with p{r) being 
the nuclear matter density. The parameterization we adopted in the numerical computations 
as well as several other technical details is given in Appendix [Bl 

The cross section can be re-expressed in terms of the gauge boson self-energy in nuclear 
medium, it is readily done by a substitution: 

L,^W^''{p{f)) = --^ [Vn'^^(g,p(r-))] . (3.2) 

TT 

The polarization tensor 11^^ has a dimension of (energy)^. After multiplying it by an appro- 
priate external couplings and performing the spatial d^r integration one gets a representation 
of the gauge boson self-energy. It can be evaluated by adding contributions from Feynman 
diagrams representing various processes, with nucleon loops having momentum cutoffs given 
by local Fermi momentum 

A dominant SPP part is in the many body language denoted as Iplhlrc (contributions 
from 2p2/il7r and more complicated final states is assumed to be small): there is one pion 
and one nucleon- hole pair (Iplh) in the final state. The corresponding contribution to 
polarization tensor can be represented as: 

ISO ^ ^ ' ^ ' 

The hadronic tensor A^pi^i^ is defined as: 

K^Hi. = J2(^'(P^'>(^)\^"c\N{p,s)){{N\p\s')n{k)\j:jN{p,s)))^ (3.4) 

s,s' 

In (13. 3p Gn denotes the nucleon propagator: 

^^^^ pO + E{p)+te\p^-E{p)-te p^-E{p)+te) ^'^ 

with UnIp) being the occupation numbers for nucleon of isospin A^. In the Fermi gas (FG) 
model nN{p) is a Heaviside step function 0(|p| — kp). The {N'{p' ,s')TT{k)\j^^\N{p,s)) 
are transition amplitudes between initial nucleon state with spin s and four-momentum p 
and final state containing pion with four-momentum k and nucleon with four-momentum 
p = p + q — k and spin s'. After inserting the nucleon propagators into polarization tensor 
we obtain the following expression: 

--^(nCiM.V) = Yl /(^ /(2vr)3 8E^{k)E{p)Eip') ^^'^^ 

[nN{p){l-nN'{p')) + nN'{p'){l-nNip))] 
6iEip') - q'^ + E^-Eip))TT[A'^;,,,^ip,q,k)] V 



with the nucleon energy E{p) = \/fP + and the final pion energy E^(k) = \Jk? + m^. 
Taking into account, that the pion may carry a charge and the nucleus atomic number can 
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be changed, one can establish the threshold corrected energy transfer (as for the quasielastic 
peak): 



In this way one accounts for the difference of rest masses of isotopes by subtracting the 
rest mass difference Qcorr. and different Fermi levels of protons and neutrons. We substitute 
qO _s. everywhere in the hadronic part of the polarisation tensor. An alternative approach 
for nuclear the binding energy is used by |JJ and shortly explained in Appendix [Bl In the 
isospin symmetric nuclei, like ^^C the exchange part of cross section given by the terms with 
n7V'(p')(l ~ ^n{p)) is negligibly small and thus we neglect it. 

3.1 Dynamics of Single Pion Production 

The dynamics is defined by a set of Feynman diagrams (Fig. [1]) with vertices determined 
by effective chiral field theory |14) . The same set of diagrams describes also pion electro- 
production, with the exception of pion pole (PP) diagram, which is purely axial. 

a) T ^ ^) T-n ^ iboson / 



Figure 1: Basic pion production diagrams from [Tl]: a) Delta pole (AP), b) crossed Delta 
pole (CAP), c) contact term (CT), d) nucleon pole (NP), e) crossed nucleon pole (CNP), 
f) pion-in-fiight (PIF), g) pion pole (PP) 

After performing summations over nucleon spins we can rewrite the hadronic tensor as: 



— <f — Qcorr + A^;^, A^;^ = Ep — Ep . 



(3.7) 




Tr + M)s^(y + M)7°s''S°] . 



(3.8) 



Reduced current matrix elements correspond to weak transition amplitudes: 



{N\p\ s')n{k) N{p, s)) = u,{p')s^u,{p). 



(3.9) 
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They are calculated to be (see: [1A\ 



sip = iC^p—cosQc ,^ , ..2_,,2 ^,,.^u^^.^2, (3-10) 



4ap = ^gcAP— cos 9c _ _ ^ , _ (3.11 



{p - kf -Ml + iMAT{{p - kf 



'^^""71^ ^ iP + ,/ ccN{<l)F.{k - q) (3.12) 

ScNP = -^C^c^vp^^ cos9cJ^cjv(9) _ _ jvf2 ^ ~ ^) (3.13) 

= -^C^CT^ COS ec7"[^7AFcT(9')7'-i^p((g-^)')]i^^(^-g) (3.14) 

M -r^ n , 2^ 2M(2/c^-g)75 ,^ 

= ~'^^'^7f^ cos9cF;,^(g ) _ _ F^(A; - g) (3.15) 

= -^Cpp^cosQcFp{{q-kf)-f^ (3.16) 

We use the convention of [14j. In our notation /* = 2.16 is the vriVA coupling constant. 
This value is slightly larger than 2.14 used in [1A\. With our choice free A(1232) width is 
0.118 [Gel^]. The values of axial couplings are standard: = 1.267 and = 93 [MeV]. We 
use averaged masses for nucleons and pions: M = |(Mn + Mp), = |(^7r+ + frij^- + m^o) 
with the values given by Particle Data Group [IB]. For the A-resonance contributions 
we assume Ma = 1.232 [GeV]. In the Delta pole (AP) and crossed Delta pole (CAP) 
amplitudes Papipii) and rA(s) denote the Rarita-Schwinger projection operator on spin-| 
states and free A — )■ t^N decay width. By T^^{p, q) we denote the A electroweak excitation 
vertex. We will give more details about the A propagator and decay width in the next 
subsection. The electroweak excitation vertex as well as the set of vector and axial form 
factors is described in Appendix [Dl For the nucleon weak currents present in fl3.12p and 
(I3.13P we use the standard vector-axial prescription: 

fccNiq) = V^{q)-A^{q) (3.17) 

V - q^ 

From the conserved vector current (CVC) hypothesis one can also get constraints on 
form factors of contact term (CT) and pion-in-fiight (PIF) diagrams: 

FpiFiQ') = Fct{Q') = F^iQ') (3.18) 

We choose the same nucleon form-factors as in [14j. Details are described in the Appendix 
[Cl Our current matrix elements contain a virtual pion form factor -F^(A; — q) coming from 
the PIF term, where the W boson interacts with a virtual pion with momentum a = k — q. 
The CVC forces one to include it in several other background terms. is assumed to have 
a monopole form: 

F.{a) = j^,_J ; K = 1.25[GeV]. (3.19) 
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The p- meson form factor Fp{a) = jz^rj^'^ = 0.7758 [GeV]. has been introduced in the 
PP term by the authors of [14j in order to account for the p-meson dominance of ttttNN 
couphng. Because of the partially conserved axial current (PCAC) hypothesis it has been 
also introduced in the axial part of CT. For each physical pion production channel there is a 
set of isospin Clebsch-Gordan coefficients Ci. They are listed in Tab. [H In the carbon cross 



Table 1: Charged current isospin coefficients of f l3.10tl3TT6|) 



Process 


AP 


CAP 


NP 


CNP 


CT, PIF, PP 


i^i + p ^ + + p 




073 





1 


1 


Ui + n ^ + 71^ + p 


-v/2/3 


0/3 


Vl/2 


-Vl/2 


-V2 


z/; + n — 7- + 7r+ + n 


0/3 


x/3 


1 





-1 


17/ + n — 7- /+ + TT^ + n 




0/3 





1 


1 


17/ + p + 7r° + n 


v/2/3 


-0/3 


-0/2 


v/1/2 


V2 




0/3 




1 





-1 



section computations we sum up contributions from protons and neutrons in the incoherent 
way. 



3.1.1 A(1232) Decay Width and Propagator 

The ttNA interaction is decribed by the Lagrangian: 

C^NA = —'^S\d^$)^ + h.c. (3.20) 
This results in the following formula for free vacuum A — vrA^ decay width: 

p.ac. _ 1 kUEN,cm + M) 

iA^iv. - ^2vrm2 W ^^'^^^ 

It is worthy to notice, that the authors of [1] and use: 

1 r'kum 

In the above formulae cm denotes the A center of mass frame. 
The default A propagator is given by: 

< (PA, A/a) = -(/a + M^) (^g" - l^^^ - + ^.Plll^^] . (3.24) 

In the above equation is the projection operator on spin-| states with p^ being the A 
resonance 4-momentum and the free resonance decay width given in (I3.22p . 
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3.1.2 A Self-Energy 

The A(1232) isobar exhibits a strongly medium-dependent behavior due to the possibility 
to decay into a pion-nucleon pair. The free resonance decay width gets decreased because 
of Pauli blocking. Assuming a uniform distribution of decay pions in the A rest frame, the 
Pauli blocking factor is calculated to be: 

+ \p^.\kcm- EpW /o orN 

^(Pa> \VM^Ef) = r^—r (3-25) 

\PA I Kcm 

and 

^ f = f{pI, |pA|,^F)r'"^^(s). (3.26) 

But inside nucleus other decay channels are opened: the two- and three-nucleon absorp- 
tion. The net effect is an overall increase of the A width. 

rr-(s) ^ 2(ifA + iET'^n = f A - 2(53SipiM. + '^^2p2h + Q^3p3h - ^3?Sa) (3.27) 

In Oset parameterized this width as a functions of either the incoming pion kinetic 
energy x = ^ or the real photon energy and the local density of nuclear matter. We use 
his approach in our computations. It is necessary to translate the Oset results obtained in 
the kinematical situations of real photon or pion scattering to the situation of virtual boson 
interaction. It was assumed that the Oset functions: 

- ^5Ea = CipihiAp/PoT + C2p2hip/pof + CsMp/PoV (3.28) 

(all Cx and a, P, 7 are the functions of photon energy or pion kinetic energy) are in a good 
approximation the functions of the average invariant hadronic system mass. The relations: 

together with W"^ = {pn + qY allow us to translate the virtual boson into one of the available 
parameterizations. For the real part of self-energy we use the same prescription as in |17] : 

3?(S) =40^[MeV] (3.30) 

This prescription neglects different renormalizations of the longitudinal and transverse A 
response functions in the nuclear medium, but for our purpose it is sufficient. 

The main problem in using these prescriptions in model of [14j comes from the fact 
that Sa is calculated using nonperturbative effects not included in tree-level diagrams of 
(13. 10113. 16|) . All of them contain simple single pion interaction vertex. Thus we modify only 
the widths in denominators of AP diagram by substituting: 

^ F (3.31) 



Pi - Mi + ^MAr—(.) pl-Ml + zMa 



r - 2(5>Sa - «3ftS 



The many-body correction to the SPP through A resonance ^JSipiftiTr and cross sections 
for multinucleon channels connected to 9=S2p2/i and 53X3^3/1 can be accounted for by changing 
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the |APp contribution (ISTTOj) . It can be done by substituting it by a full A resonance 
production cross section: 



dE'dVL' 



16TT^\k\. 



E{p){Ma+W) (H^_(MA+3ftSA))2+(if-53SA)2 



(3.32) 



Tr 



7°r"^%°p„f (PA)r^^(/ + M) 



The approximation comes from the nonrelativistic expansion in the A propagator p^ — + 
zMaTa ~ (Ma + W^)(PF — Ma + f Ta). As for the isospin dependence: for electrons proton 
and neutron get the same factor of 1; for neutrinos/antineutrinos protons/neutrons get a 
factor of 3 because of the Clebsch-Gordan \/3 in the weak A excitation vertex. 



4 Numerical procedures 

The full integration of cross-section within LDA (as given in Eqs 13.11 and 13. 2p even with an 
assumption of spherically symmetric nuclear matter distribution and on-shell nucleons would 
require performing six nested integrals. For a small (9(10) number of integration points in 
each of them we would need to evaluate O{10^) points in the numerical integration procedure 
to obtain just one point in the triple-differential cross-section. Thus the authors of [T] 

assumed the nucleon momentum to be an average one in local Fermi sea, = \J\p^{f)- 

Furthermore p is assumed to be orthogonal to the (g, k) plane. Within this approximation 
the number of nested integrals is reduced by 2: 

1 P d^k 1 

iso ^ ^ V 7 

/ • d^p [nN{p){l-nN>{p')) +nN'{p'){l-nN{p))] 
J (27r)3 AE{p)E{p') 

5{E{p')-q'^ + E^-E{p)) 

The integral over d^p can now be performed analytically, giving a result proportional to the 
Lindhard function. There are severe shortcomings of this approximation and we loose a lot 
of precision. One example is the threshold behaviour of the pion production cross section. 
The hadronic tensor is described by an averaged invariant pion-nucleon mass. Thus the 
physically meaningful tensor is obtained, when 

{W^) = M^ + 2(E,v)g° + > (M + m^f. (4.2) 

The above mentioned condition is important for nucleon pole (NP) diagram, for which 
an unphysical W may give rise to a sing ularity at {{p) + qf = M^. This requires an 
additional cutoff in the acceptable kinematics, which sometimes moves up the threshold for 
pion production process in an artificial way. 

However, the six dimensional integration can be performed using Monte Carlo tech- 
niques. There exist several available algorithms for that. We have chosen the Vegas algo- 
rithm implemented in GNU Scientific Library (GSL) for C/C++ compilers ^U\. It is efficient 
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Figure 2: Difference between the exact cross section calculation from this paper and ap- 
proximations used in 

enough to perform 8-dimensional total cross section integration in a reasonable time using 
only O{10^) points. This solves the threshold problem caused by averaged hadronic tensor 
with averaged W"^. 

In order to show the difference between the exact calculation and the approximation 
adopted in ( 14. ip we calculated a sample double-differential electron neutrino cross section 
off carbon. The results are shown in Fig. [2] for neutrinos (top) and for antineutrinos 
(bottom). The curves calculated using (14. ip are quite different from those calculated without 
approximations. 

For total cross-section both approaches: exact and approximate give similar results, as 
one can see in the Fig. |3l In the case of antineutrino charged pion production there is 
a systematic difference between our calculation and approximated results, but it is rather 
small. Thus we find the approximation fl4.ip sufficient on the level of total cross-sections. 
However, in what follows we will always use the exact calculations. 
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Figure 3: Difference between exact cross section calculation from this paper and approxi- 
mations used in [Ij. 

5 Results 

5.1 Importance of background terms 

Fig. m shows importance of background terms for pion production on a set of 6 free protons 
and 6 free neutrons. The curves describes ratios of cross sections coming from only Delta 
pole diagram to the cross section calculated with all the background diagrams (and their 
interference terms) included in computations. 

We see that especially for the lower neutrino energies, below 500 MeV, the background 
contribution is very important. The background terms are more relevant for antineutrinos 
than for neutrinos and for the vr^ production than for a charged pion production. 

5.2 Importance of in-medium effects 

Fig. [5] shows an impact of the in-medium effects on the pion production. We plotted a 
relative modification of the free nucleon cross section (six free protons and neutrons but 
with the background contribution included) caused by the in-medium effects. In almost 
all of the cases the in-medium effects leads to a significant decrease of the total cross- 
section. For the electron (anti)-neutrinos 400 MeV is far from the SPP reaction threshold 
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Figure 4: Ratios of the total cross sections for z/^ and z/^ SPP reactions on carbon calculated 
with a model without the background terms to the full model of this paper. 

and in those cases we see an almost constant reduction of the cross section on the level 
of 30 — 40%. There is an interesting difference in shapes between electron neutrinos and 
antineutrinos, see Fig. [5l The latter exhibits a smooth drop of in-medium reduction with 
growing neutrino energy. In the case of muon neutrinos and antineutrinos near the pion 
production threshold {E,^ < 0.5 [GeV]) the cross section is less affected by nuclear effects. 
For TT^ production channel and = OA [GeV] it even seems to be slightly enhanced. This 
happens due to nucleon Fermi motion which dominates other effects in that kinematical 
region. This is not the case for 7r° production by antineutrinos. There exists a correlation 
between the nonresonant background contribution and the cross section reductiondue to in- 
medium effects. Shapes of the reduction ratios in neutrino vr'^ and antineutrino 7r~ channels 
are almost the same, so are background contribution shown in Fig. |H In general, the 
more cross section comes from background and interference terms, the smaller is the near 
threshold effect. For the larger muon neutrino/antineutrino energies E > 0.6 [GeV] we see 
again an almost uniform reduction of the cross section of the order of 30%. 

5.3 Total cross sections 

We compared predictions from our model with the recent MiniBooNE pion production 
data. MiniBooNE, unlike K2K, published their results in a form of absolutely normalized 
cross section and not as a ratio to CC inclusive cross sections. We performed calculations 
with our model of the total cross-sections on CH2. A direct comparison with the data 
is not straightforward because MiniBooNE reported the cross sections for pions in the 
final state after leaving nucleus (in a case of neutrino-carbon scattering) with all the FSI 
effects included. The pion FSI effects can be evaluated within a cascade models like those 
implemented in Monte Carlo event generators. Our model is not yet an ingredient of any 
MC generator and we tried to estimate an impact of FSI effects using the results of MC 
comparison study published in [2l]. We approximate the relevant probabilities as: 



P(7r0 ^ 7r°) = 67%, P(7r° ^ 7r+) = 5% 
P(7r+ ^ 7r+) = 69%, P(7r+ ^ 7r°) = 5% 



(5.1) 
(5.2) 
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Figure 5: Impact of nuclear effects on SPP off Plots show {crfree~ (^medium) / c^free x 100% 



The results for the cross section with and without FSI are plotted in Fig. [6l 

In the case of charged pion production we obtained a quite good agreement with the data 
up to the neutrino energy of around 0.8 [GeV]. In the case of charged- current l7r° production 
both free and in-medium cross sections with our model are too small, and the discrepancy 
becomes larger with increasing neutrino energy. FSI introduce large modifications for the 
Tr"*" channel. In the vr'' channel an effect of absorption of 7r° is partially compensated by 
a fraction of initial 7r+ events, that end up as 7r° due to charge exchange reaction inside 
nucleus. It is important to point out that in the case of CCvr^ reaction also the computations 
of other theoretical groups give results well below the measured cross section j22) . 



5.4 Ratios of muon to electron (anti-) neutrino cross sections 

In neutrino oscillation appearance experiments it is very important to calculate precisely 
ratios of muon and electron neutrino cross sections. Even in a presence of a near detector 
and with full understanding of initial muon neutrino flux a good knowledge of the ratios (and 
their dependence on neutrino energy) is crucial for a correct identification of the oscillation 
signal. 

In Fig. [7] we see that the ratios calculated with the complete model are slowly increasing 
functions of the neutrino energy. In the case of antineutrinos there is a small difference 
between tt~ and vr'' production: in the first case the ratio is slightly lower. On the contrary, 
we obtain almost the same ratios both for 7r+ and 7r° production by neutrinos. 

It is important to know how well the ratios are calculated when much simpler models 
are used, which is a case in MC event generators. 
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Figure 6: Charged and neutral pion production cross sections on CH2 for full model of this 
paper plotted against the data from by [9] and [10]. 



100 
80 

40 
20 




T ' 1 ' r- 



1 ' 1 ' r 



neutrinos 

71^ 

71° 



*♦ ** 



antineutrinos 
n 

71° 



0.4 0.5 0.6 0.7 0.8 0.9 1 0.4 0.5 0.6 0.7 0.8 0.9 1 
Ev [GeV] E, [GeV] 



Figure 7: Ratios of muon to electron (anti) neutrino total SPP cross sections on ^^C for the 
full model of this paper. 



Fig. [8] shows an impact of the background terms on the 7r° production ratios. We 
compared two situations: the full model and the model without background contributions. 
We see that the results are significantly different only in the case of antineutrinos. For 
lower neutrino energies one obtains much smaller ratios while using pure resonant SPP 
mechanism. For neutrinos these differences are negligible. 

Fig [9] show an impact of A self-energy on the ratios. We compared two situations: the 
full model and the model without A self-energy. We see that the negligence of the A self- 
energy has almost no impact on the considered observable. We conclude, that in order to 
describe well the anti-muon to anti-electron neutrino cross section ratio it is important to 
include the nonresonant background, but not necessarly the A self-energy. 



5.5 Pionless A decays 

An interesting feature of the model we discuss is that we obtain a contribution to the cross 
section coming from pionless A decays. This is a part of the meson exchange current (MEC) 
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Figure 8: Ratios of muon to electron (anti-)neutrino total SPP cross sections on ^^C for the 
full model of this paper and for the resonant SPP only (without the background terms). 
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cross section which has recently attracted a lot of attention [23] • There is a lot of evidence 
that MEC mechanism is responsible for a large CCQE axial mass measurement reported 
by the MiniBooNE collaboration. Theoretical microscopic computations always include 
pionless A decays as a part of the calculated effect. Some MC event generators (NEUT, 
NUANCE) assume a constant fraction of the pionless A decays and we find it interesting 
to check how well this assumption is satisfied in our model. 

The fractions of the pionless decays and their dependence on the neutrino energy and 
species are shown in Fig. [TOl There is no difference between neutrinos and antineutrinos, 
because we include only the np — nh mechanism coming from the resonant diagrams. The 
fraction of pionless A decays is very large for the energies below 500 MeV. For the larger 
energies it exhibits a smooth energy dependence, dropping down to 20% sX E^^ = 1 [GeV]. 
It is clear that for experiments with a large fraction of neutrinos with energies below 1 [GeV] 
one can not consider the investigated quantity to be constant. 

The total pionless A decay cross section may be treated as a lower bound for the np — nh 
contribution. One has to keep in mind, that there are many more sources of np — nh final 
states, which can be built from diagrams (13. 10113. 16|) but are not considered in this paper. 
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Figure 10: Fraction of the pionless A decays to the resonant SPP production cross section 
{o-pioniess a)/(^spp res. X 100% iu ^"^C for z/g and i/^. 

The total cross section coming from A ^ np — nh decays can be seen in Fig. [11] The 
charge current quasielastic (CCQE) contribution has been calculated with NuWro neutrino 
event generator [21] with Ma = 1-05 [GeV]. The np — nh contribution coming from pionless 
A decays is small compared to CCQE and SPP dynamics (around 10-15% of the first), but 
it is non-negligible. Alone, it cannot explain MiniBooNE's large axial mass measurement 
because that this is only a part of the np — nh cross section. 

6 Discussion 

It is not east to understand why the A self-energy leads to the substantial reduction of 
the cross sections. As explained in the introduction the in-medium A spectral function 
was included only in the AP diagram. It is difficult to conclude, whether the cross section 
reduction is a genuine physical effect or rather an artifact of the approximate treatment of 
background terms and in-medium A self-energy. The pure background contribution (36 out 
of 49 combinations from f l3.10p -( l3.16p ) is not affected by the presence of nuclear matter. In 
the A-background interference terms (12 combinations) the in- medium effects enter only 
through AP diagram, thus are included only partially. A complete in-medium dressing is 
present only in the pure A contribution, reducing its size significantly. A verification of 
the model prediction can come only from the evaluation of the non-perturbative in-medium 
effects for all the genuine amplitudes (28 independent terms) which is a very difficult task 
to achieve. 

We arrived at a reasonable agreement with the MiniBooNE CCvr"*" production but our 
model underestimates CCvr^ cross section. There can be several reasons for that. The 
first one can be approximations discussed in the previous paragraph. For neutrino energies 
around 1 GeV one should include also contributions from heavier resonances absent in 
our computations. It is also possible that the lTT2p2h process contributes with a larger 
cross section than it is generally expected. An apparent excess of the CCItt^ cross section 
with respect to several theoretical models predictions is an interesting topic for the further 
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Figure 11: Total charge current cross sections on ^^C for: quasielastic scattering, SPP, and 
pionless A decay. 

research. 

We investigated a possible impact on predictions from the model coming from different 
descriptions of the A resonance width. For example the authors of [25] use: 

Tm-s{W) = 118[Me\/] ■ ^^^^^^ (6.1) 

The term ^^"p ^2 is a so-called Blatt-Weisskopf centrifugal barrier. In this manner one 
accounts for t£e phenomenological knowledge about decay txN system angular momentum, 
which is absent in the Lagrangian f l3.20p . Furthermore, A(1232) is not a stable particle. 
One can account partially for the off-shell effects by replacing the propagator in DP term 
by 

r.,(^. = ^3"/2(pa) _ (XA + iy) .... 

^^^^ Ml + iWT^{pl) pl-Ml + iWT^ipl)'' ^^-"^ 




3 ' ' 3 3 W 

This convention is used by |26] together with the Manley-Saleski decay width. Thus while 
we use (16. ip we also replace (I3.23P by (16. 2p . In order to stay consistent, after changing the 

width (I3.2ip with (16.11) in (16. 2p one has to multiply the whole expression by \J ^rA{^ ^' 
will compensate for the fact, that our current has a decay vertex defined by (I3.20p in the 
numerator, which leads to the width (13.210 . 

On the level of total cross sections we find the difference between two AttN decay 
descriptions negligible. This is illustrated in the Fig. [T2] where we plot again the muon to 
electron (anti-)neutrino total vr^ production cross section ratios and we compare the default 
and the Manley-Saleski A description. Both descriptions lead to almost identical results. 
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Figure 12: Ratios of muon to electron (anti) neutrino total SPP cross sections on ^^C for 
the full model of this paper with A width described by (13.211) and (16. ip . 
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Finally, we would like to address the question: how much does the numerical approxi- 
mation (14. ip affect muon to electron neutrino cross section rates. We have already shown, 
that the exact integration does not change much total cross sections. This is illustrated in 
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Figure 13: Ratios of i/^ to z/g neutrino CCvr" SPP cross sections on ^^C calculated with the 
full model of this paper and with approximations used in [H]. 

Fig. [131 where we have plotted V^/Ve ln^CC cross section ratios. Differences are only for 
energies < 550 [MeV] and at Eu = 500 [MeV] it is about 4.3%. 
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A Notation and conventions 

We adapt the conventions from Bjorken-Drell {g'^'^ = (+, — , — , — ) etc.), the only difference 
is in the Dirac spinor normahzation: 

^Us{p)us{p) = {/ + m) (A.l) 

s 

which is convenient for our calculations. 

B Nuclear Matter Density Parameterization 

We took the harmonic oscillator density profiles from |27| : 

p(f) = po (1 + « {r/Rf) exp [- [r/Rf] (B.l) 

with corrections to parameters a and R calculated in [28J. These parameters are slightly 
different for protons and neutrons. The local Fermi momentum is calculated from relation 
kp{r) = (37r^p(r*)^) 3 . Authors of [Tj subtract Fermi kinetic energy from nucleons inside 
medium E{p) — > E{p) —Tp. In this manner they account for the binding effects. 

C Nucleon Form Factors 

The isospin symmetry relates the vector form factors to the electromagnetic ones: 

Frm=F!m-F-m- (c.i) 

For the electromagnetic form factors we use the parameterization of Galster et al. [29]: 



with /ip = 2.792.847, /i„ = 1.913043, A„ = 5.6, r = ^ and = 0.843 [GeV]. We assume 
the axial nucleon form factor in a dipole form: 

with = 1.267. 
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D A(1232) Form Factors 

The most general electroweak A excitation vertex is given by: 
r"^(p,g) = 



Tra/1 _ ACtfJ. 

^3/2 ^3/2 

a 



(D.i: 



CY 



M 



M2 



+ ^{g'^'^q-P-qy^+g'^^c^ 



+ 



^{g'''M'-qy') + %{g'''q-{p+q)-q''{p+qm 



Ci 



M2 



The and Cf" vector and axial form factors determine WN [\ transition. For the vector 
form-factor set we use the parameterization of |30| : 

2.13 1 



cl 



(D.2) 



-1.51 



1 



:i + Q7M2)2 1 + QV4M2 
0.48 1 



;i + QVM2)2 1 + Q2/0.776M2 



with My = 0.84 [GeV]. The CVC implies that CX = 0. The axial part is dominated by 
contribution. We use a dipole approximation, in which 

cm 



(D.3) 



We use a default value of the A axial mass M^a = 1-05 [GeV]. The default value of C^(0) 
is obtained from the Goldberger-Treiman relations [3lj : 



V 3m^ 



(D.4) 



which is somewhat higher, than what is used in \\^: C^ifd) ~ 1.15. The authors of [30 

-A/^2N _ CiiO) 1 



and P] use C^{Q'^ 

transition, which do not allow to extract any beyond-dipole behavior, we use the simple 



f^^n 2l^fi \2^,n2l■}^fi ■ Because of big Uncertainties in axial iVA 



dipole form (lD.3p . We include Adler [32] relation for C4 , i. e.: 



Furthermore, from PCAC hypothesis one can determine: 

M2 



ictm- 



ml + 

The form factor is considered to be negligibly small, thus we set C^iQ"^) = 



(D.5) 



(D.6) 
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E Total Cross Section Tables 

In the following section we present our results in form of the tables. 
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Table 2: Total cross-sections in the ^^C(z/e,e ) scattering in 10 [cm^]. 
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1.0477 
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Table 3: Total cross-sections in the ^^C(z/e, e"*") scattering in 10 [cm^]. 
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Table 4: Total cross-sections in the ^^C(z/^,/i ) scattering in 10 [cm^]. 
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0.3425 


0.7725 


0.1545 


1.0116 


0.2244 


0.3667 


0.65 


1.6128 


0.3226 


1.9898 


0.4655 


1.7418 


0.3484 


2.0115 


0.4358 


1.0382 


0.2076 


1.3167 


0.2904 


0.4417 


0.70 


2.0407 


0.4081 


2.4533 


0.5704 


2.1309 


0.4262 


2.4195 


0.5230 


1.3044 


0.2609 


1.6114 


0.3554 


0.5099 


0.75 


2.4628 


0.4926 


2.8998 


0.6718 


2.4985 


0.4997 


2.7860 


0.6017 


1.5630 


0.3126 


1.8779 


0.4188 


0.5705 


0.80 


2.8596 


0.5719 


3.3105 


0.7661 


2.8390 


0.5678 


3.1156 


0.6740 


1.8092 


0.3618 


2.1372 


0.4754 


0.6238 


0.85 


3.2263 


0.6453 


3.6855 


0.8533 


3.1519 


0.6304 


3.4359 


0.7380 


2.0404 


0.4081 


2.3644 


0.5294 


0.6700 


0.90 


3.5638 


0.7128 


4.0288 


0.9345 


3.4382 


0.6876 


3.6895 


0.8048 


2.2549 


0.4510 


2.5930 


0.5776 


0.7098 


0.95 


3.8729 


0.7746 


4.3430 


1.0100 


3.6954 


0.7391 


3.9564 


0.8624 


2.4532 


0.4906 


2.7903 


0.6250 


0.7439 


1.00 


4.1550 


0.8310 


4.6313 


1.0805 


3.9245 


0.7849 


4.1873 


0.9113 


2.6352 


0.5270 


2.9678 


0.6661 


0.7730 



Table 5: Total cross-sections in the ^^C(z/^, /i"*") scattering in 10 ^® [cm^]. 
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0.1649 
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0.1817 
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0.0748 
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0.1765 


0.0568 
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0.65 
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0.3546 
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0.2432 
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0.70 
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0.0662 


0.4531 


0.1286 


0.1981 


0.0396 


0.3167 
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0.1244 
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0.0633 


0.4871 


0.1509 
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0.85 


0.5841 


0.1168 


0.8424 


0.2883 


0.5962 


0.1192 


0.7937 


0.2237 


0.3823 


0.0765 


0.5811 


0.1783 
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0.90 


0.6792 
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0.9733 


0.3321 
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0.8926 


0.2685 





1709 



References 

[1] J. Nieves, I. Ruiz Simo, and M. J. Vicente Vacas. Phys. Rev. C, 83:045501, 2011 

[2] A. A. Aguilar-Arevalo et al. [MiniBooNE Collaboration]. Phys. Rev. D, 81:092005, 
2010 

[3] K.M. Graczyk and J.T. Sobczyk. Phys. Rev. D, 77:053001, 2008. erratum-ibid. D79 
(2009) 079903 

[4] M. Martini, M. Ericson, and G. Chanfray. Phys. Rev. D, 85:093012, 2012. ; 

J. Nieves, F. Sanchez, I. Ruiz Simo, and M. J. Vicente Vacas. Phys. Rev. D, 85:113008, 

2012 ; 

O. Lalakulich and U. Mosel. arXiv: 1208.3678 [nucl-th] 

[5] L. Tiator, D. Drechsel, S. S. Kamalov, and M. Vanderhaeghen. Eur. Phys. J. ST, 
198:141, 2011 

[6] D. Drechsel, S. S. Kamalov, and L. Tiator. Eur. Phys. J. A, 34:69, 2007 

[7] A. Rodriguez and L et al. [K2K Collaboration] Whitehead. Phys. Rev. D, 78:032003, 
2008 

[8] C. Mariani et al. [K2K Collaboration]. Phys. Rev. D, 83:054023, 2011 

[9] A. A. Aguilar-Arevalo et al. [MiniBooNE Collaboration]. Phys. Rev. D, 83:052007, 
2011 

[10] A. A. Aguilar-Arevalo et al. [MiniBooNE Collaboration]. Phys. Rev. D, 83:052009, 

2011 

[11] S. J. Barish, M. Derrick, T. Dombeck, L. G. Hyman, K. Jaeger, B. Musgrave, 
P. Schreiner, and R. et al. Singer. Phys. Rev. D, 19:2521, 1979 

[12] G. M. Radecky, V. E. Barnes, D. D. Carmony, A. F. Garfinkel, M. Derrick, E. Fernan- 
dez, L. Hyman, and G. et al. Levman. Phys. Rev. D, 26:3297, 1982. [Erratum-ibid. D 
26 (1982) 3297] 

[13] T. Kitagaki, H. Yuta, S. Tanaka, A. Yamaguchi, K. Abe, K. Hasegawa, K. Tamai, and 
H. et al. Sagawa. Phys. Rev. D, 42:1331, 1990 

[14] E. Hernandez, J. Nieves, and M. Valverde. Phys. Rev. D, 76:033005, 2007 

[15] O. Lalakulich, T. Leitner, O. Buss, and U. Mosel. Phys. Rev. D, 82:093001, 2010 

[16] E. Oset and L. L. Salcedo. Nucl. Phys. A, 468:631, 1987 

[17] S. K. Singh, M. J. Vicente- Vacas, and E. Oset. Phys. Lett. B, 416:23, 1998. [Erratum- 
ibid. B 423 (1998) 428]; 

S. Ahmad, M. Sajjad Athar and S. K. Singh. Phys. Rev. D, 74:073008, 2006; 
M. Sajjad Athar, S. Ahmad, and S. K. Singh. Nucl. Phys. A, 782:179, 2007; 
S. K. Singh, M. Sajjad Athar, and S. Ahmed. AIP Conf. Proc, 967:182, 2007; 
M. Sajjad Athar, S. Chauhan, and S. K. Singh. J. Phys. G, 37:015005, 2010; 
M. Sajjad Athar, S. Chauhan, and S. K. Singh. Eur. Phys. J. A, 43:209, 2010 



26 



[18] J. Beringer et al. [Particle Data Group Collaboration]. Phys. Rev. D, 86:010001, 2012 
[19] A. Gil, J. Nieves, and E. Oset. Nucl. Phys. A, 627:543, 1997 

[20] M. Galassi et al. GNU Scientific Library Reference Manual (3rd Ed.). Network Theory 
Ltd., 2009 

[21] M. Antonello, V. Caracciolo, G. Christodoulou, J. Dobson, E. Prank, T. Golan, V. Lee, 
and S. et al. Mania. Acta Phys. Polon. B, 40:2519, 2009 

[22] O. Lalakulich, K. Gallmeister, T. Leitner, and U. Mosel. AIP Conf.Proc, 1405:127, 
2011 

[23] J. Marteau, J. Delorme, and M. Ericson. Nucl. lustrum. Meth. A, 451:76, 2000. ; 
M. Martini, M. Ericson, G. Chanfray, and J. Marteau. Phys. Rev. C, 80:065501, 2009; 
J. Nieves, L Ruiz Simo, and M. J. Vicente Vacas. Phys. Lett. B, 707:72, 2012 

[24] C. Juszczak, J. A. Nowak, and J. T. Sobczyk. Nucl. Phys. Proc. SuppL, 159:211, 2006. 
J. A. Nowak. Phys. Scripta T, 127:70, 2006; 
C. Juszczak. Acta Phys. Polon. B, 40:2507, 2009 

[25] D. M. Manley and E. M. Saleski. Phys. Rev. D, 45:4002, 1992 

[26] T. Leitner. Neutrino-Nucleus Interactions in a Coupled- Channel Hadronic Transport 
Model. PhD thesis, Justus-Liebig-Universitaet Giessen, Fachbereich 07 (Mathematik 
und Informatik, Physik, Geographic), Institut fuer Theoretische Physik, 2007 

[27] H. De Vries, C. W. De Jager, and C. De Vries. Atom. Data Nucl. Data Tabl., 36:495, 
1987 

[28] C. Garcia-Recio, J. Nieves, and E. Oset. Nucl. Phys. A, 547:473, 1992 

[29] S. Galster, H. Klein, J. Moritz, K. H. Schmidt, D. Wegener, and J. Bleckwenn. Nucl. 
Phys. B, 32:221, 1971 

[30] O. Lalakulich, E. A. Paschos, and G. Piranishvih. Phys. Rev. D, page 014009, 2006 
[31] M. L. Goldberger and S. B. Treiman. Phys. Rev., 110:1178, 1958 
[32] S. L. Adler. Annals Phys., 50:189, 1968 



27 



100 
80 
60 
40 
20 




/ a(v^)/cy(Ve) 



.+ 



71^ exact 



71 approx. 



0.4 0.5 0.6 0.7 0.8 0.9 

Ev [GeV] 



